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Abstract — This paper proposes a new method for optimizing 
frequency-hopping ad hoc networks in the presence of Rayleigh 
fading. It is assumed that the system uses a capacity-approaching 
code (e.g., turbo or LDPC) and noncoherent binary continuous- 
phase frequency-shift keying (CPFSK) modulation. By using 
transmission capacity as the performance metric, the number 
of hopping channels, CPFSK modulation index, and code rate 
are jointly optimized. Mobiles in the network are assumed to be 
uniformly located within a finite area. Closed-form expressions 
for outage probability are given for a network characterized 
by a physical interference channel. The outage probability is 
first found conditioned on the locations of the mobiles, and 
then averaged over the spatial distribution of the mobiles. 
The transmission capacity, which is a measure of the spatial 
spectral efficiency, is obtained from the outage probability. The 
transmission capacity is modified to account for the constraints 
of the CPFSK modulation and capacity-approaching coding. 
Two optimization methods are proposed for maximizing the 
transmission capacity. The first is a brute-force method and 
the second is a gradient-search algorithm. The results obtained 
from the optimization shed new insight into the fundamental 
tradeoffs among the number of frequency-hopping channels, the 
modulation index, and the rate of the error-correcting code. 

I. Introduction 

Ad hoc networks comprise mobiles that communicate with- 
out centralized control or a pre-existing infrastructure. The pre- 
ferred channel access for ad hoc networks is direct-sequence or 
frequency-hopping (FH) spread spectrum. This paper focuses 
specifically on frequency-hopping spread spectrum ad hoc 
networks. Such networks are characterized by independent, 
identical, FH radios that share the same carriers and frequency 
channels, and are nearly stationary in location over a single 
hop duration. 

The first part of this paper is concerned with the analysis 
of the outage probability of FH networks, where outage 
probability is the probability that the signal-to-noise-and- 
interference ratio (SINR) falls below a predetermined thresh- 
old. By limiting the fading to be slow Rayleigh fading and 
excluding shadowing, the paper first presents an exact closed- 
form expression for the outage probability conditioned on the 
locations of the interferers. 

The interferers are assumed to be uniformly distributed in 
an annular area, where the inner radius is a minimum interferer 

This work was supported in part by the National Science Foundation under 
Award No. CNS-0750821 and by the United States Army Research Laboratory 
under Contract W91 1NF-10-0109. 



distance that could be imposed by an interference-avoidance 
protocol HI, such as carrier-sense multiple access, and the 
outer radius is the maximum distance set by the network's 
geographic footprint. By averaging over the uniform locations 
of the interferers, the spatially averaged outage probability 
is obtained in closed form. A distinguishing feature of this 
paper is that it considers networks of limited area, in contrast 
with the current popular literature, which typically assumes 
networks of infinite extent (e.g., Q, ||3l). 

The number of mobiles in the network may be either fixed 
or random. Initially, a fixed number of mobiles is assumed, 
in which case the mobile locations are a realization of a 
binomial point process (BPP). Next, it is assumed that the 
number of mobiles is Poisson distributed, in which case the 
mobile locations are a realization of a Poisson point process 
(PPP). Considering a PPP allows us to obtain results that 
are consistent with the current popular literature, and in fact, 
our results coincide with the infinite-network results of j4j, 
[5 1, |6| when we let the network boundary extend to infinity. 
However, the BPP results are of practical interest because of 
limitations to the PPP model. The most significant limitation 
to the PPP model is that it allows an unbounded number of 
users, which is not possible in a finite network. However, the 
PPP has been favored in the literature because it enables the 
use of Campbell's theorem [7], which often leads to tractable 
mathematical expressions that vastly simplify the performance 
analysis. 

Having found the spatially averaged outage probability 
under the BPP and PPP models, the paper next derives closed- 
form expressions for the transmission capacity (8), which 
is the spatial spectral efficiency; i.e., the rate of successful 
transmissions per Hz and m 2 . We propose a modification to 
the transmission capacity metric of [8| that accounts for mod- 
ulation and coding constraints. The utility of the modulation- 
constrained transmission capacity is that it can be used to 
optimize the main parameters that influence the network's 
performance. The network's performance depends on several 
parameters related to the choice of modulation and coding, and 
also depends on the number of hopping channels. It is assumed 
that the system uses noncoherent binary continuous-phase 
frequency-shift keying (CPFSK) modulation, which is the 
most common choice of modulation for FH systems [9|. The 
main parameter associated with binary CPFSK is the modula- 



tion index, which characterizes the relative separation between 
the two tones. It is furthermore assumed that the system uses 
a capacity-approaching code (e.g., turbo or LDPC), which 
allows the achievable performance to be characterized by the 
capacity of the system, under constraints of the modulation 
and noncoherent detection technique. Under the assumption of 
coded noncoherent CPFSK, the performance of the network is 
a function of three parameters: the code rate, the modulation 
index, and the number of hopping channels. 

By using the modulation-constrained transmission capacity 
as the objective function, the paper optimizes the network 
with respect to these three parameters. Initially, a brute-force 
exhaustive optimization is proposed that optimizes over a wide 
range of discretized parameters. Because the results of the 
exhaustive optimization suggest that the optimization problem 
is convex, a gradient search algorithm ifTUl is proposed that 
offers a good tradeoff between accuracy and efficiency. 

The main contributions of this paper are (1) the closed-form 
expressions for conditional outage probability and spatially 
averaged outage probability in the presence of finite-area 
networks with mobiles drawn from both a BPP and a PPP, 
(2) the development of modulation constrained transmission 
capacity as a performance metric, and (3) a method for 
optimizing the parameters associated with the ad hoc network. 
The methodology presented in this paper presents a new 
approach to the analysis and optimization of finite ad hoc 
networks and presents fresh insight into the tradeoffs among 
the number of frequency-hopping channels, the modulation 
index, and the code rate in a frequency-hopping network. 

II. Network Model 

The network comprises M + 2 mobiles that include a ref- 
erence receiver, a reference transmitter Xq, and M interfering 
transmitters Xi,...,Xm- The coordinate system is selected 
such that the receiving mobile Xq is at the origin. The variable 
Xi represents both the i th mobile and its location, and \ \Xi\\ 
is the distance from Xi to the receiving mobile. While the 
interferers can be located in any arbitrary region, we assume 
they are located in an annular region with inner radius r ex and 
outer radius r net . A nonzero r ex may be used to model the 
effects of interference-avoidance protocols (TJ. In particular, 
a nonzero r ex models an exclusion zone placed around the 
receiver, which can be realized by having the receiver send 
a short clear-to-send (CTS) packet in response to a request- 
to-send (RTS) packet sent by the transmitter. Under a carrier- 
sense multiple-access (CSMA) protocol, mobiles within dis- 
tance r ex from the receiver that overhear the CTS will suppress 
their transmission. 

Xi transmits a signal whose average received power in the 
absence of fading is Pj at a reference distance ro- At the 
receiving mobile, Xi's power is 



P> 



Piffi/(||*i||) 



(1) 



where <?j is the power gain due to fading, and /(||Xj||) is a 
path-loss function. Each gi = af, where a% is Rayleigh and 



§i unit-mean exponential, i.e. ~ £(1)- For r > r , the 
path-loss function is expressed as the attenuation power law: 



f(r) = 



(2) 



where a > 2 is the attenuation power-law exponent and r is 
sufficiently large so that the signals are in the far field. 

Channel access is through a synchronous frequency-hopping 
protocol. The hopping is slow, with multiple symbols per hop, 
which is a more suitable strategy for ad hoc networks than fast 
hopping ifTTI . An overall frequency band of B Hz is divided 
into L frequency channels, each of bandwidth B/L Hz. The 
transmitters independently select their transmit frequencies 
with equal probability. Let pi denote the probability that 
interferer Xi selects the same frequency as the source. Let 
di < 1 be the duty factor of the interferer. It follows that 
Pi — di/L and that using a duty factor less than unity is 
equivalent to hopping over more than L frequencies IfTTI . 
Assuming that di — d for all interferers, L' = L/d denotes 
the equivalent number of frequency channels. It is assumed 
that the {g.i} remain fixed for the duration of a hop, but vary 
independently from hop to hop (block fading). While the {gi} 
are independent from user to user, they are not necessarily 
identically distributed. 

The instantaneous SINR at the receiving mobile is 



7 



Po 



(3) 



Af + hPi 



where Af is the noise power and I; is a variable that indicates 
the presence and type of interference (i.e. co-channel interfer- 
ence or adjacent-channel interference). When adjacent-channel 
interference IfTTI is neglected, Ii = 1 when Xi selects the same 
frequency as Xq, and p = otherwise. It follows that Ii is 
Bernoulli with probability P[Ii = 1] = p t . 
Substituting ([]} and Q into the SINR is 



M 

E 

1=1 



(4) 



where T = r^P^/Af is the signal-to-noise ratio (SNR) 
when the transmitter is at unit distance and fading is absent, 
f2< = (-Po/P)||^i|| Q i s me inverse normalized power of Xi 
at the receiver, and VLq = ||Xo|| a . Without loss of generality^ 
we assume that ||Xo|| = 1 for the remainder of this paper. 
Furthermore, all examples and numerical results in this paper 
assume that all mobiles transmit with the same power, i.e. 
P = P for all i. 

III. Conditional Outage Probability 

Let /3 denote the minimum SINR required for reliable 
reception and J~2 = {f2 , ^a/} represent the set of inverse 

'Changing ||Xn|| is equivalent to a scaling of r ex and r ne t. 



normalized powers. An outage occurs when the SINR falls 
below (3. Conditioning on f2, the outage probability is 



en = P [7 </?|«] 



(5) 



Because it is conditioned on the outage probability depends 
on the particular network geometry, which has dynamics over 
timescales that are much slower than the fading. By defining 
a variable M 

z M = /r V^o 1 - ^ W 1 (6) 



the conditional outage probability may be expressed as 

ea = P [Z M < T- 1 ^] = F Zm (T-'ln) (7) 

which is the cumulative distribution function (cdf) of Zm 
conditioned on fl and evaluated at T" 1 . 

By defining S = f3~ 1 flQ 1 g a and Yj = Ij^jf^ -1 , j6j may be 
rewritten as M 

(8) 



Z M = S-^Y, 

<=i 



where S ~ £(/3fio)- The cdf of S is 



Fs(y) 



(9) 



where u(y) is the unit-step function. Taking into account the 
Rayleigh fading and Bernoulli {ij}, the pdf of Yj is 



where 5(y) is the Dirac delta function. 

First consider the single-interferer case (M 
of Zi is expressed as 

/•OO 

F Zl (z|n) = / F s (z + y)f Yl (y)dy 
Jo 

Substituting ^ and ( [10} into ( fTT) yields 



(10) 
1). The cdf 
(11) 



F Zl (z\n) 



-/3Q (z+y) 



= 1 - e 



-l3a z 



(1 



hi(y)dy 



(12) 



for z > 0. Using the fact that Z M - 
iteratively 

M 

F TM {z\n) = l-e-^'H 



-M-l 



Yjj, and working 



(1 -Pi)jgn H 



(13) 



for z > 0. The outage probability is found by substituting 
z = r _1 into the above expression. 

Example #1: Consider a specific network topology with 
the transmitting mobile placed one unit North of the receiving 
mobile, and fifty interferers arbitrarily placed in an annular 
region with outer radius r net = 2 and inner radius r ex = 0.25. 
The resulting network is shown in the inset of Fig. [T] The ft 
was determined by assuming a path-loss exponent a = 3 and 
a common transmit power Pj = Po- The equivalent number 
of frequency channels was set to V = 200. Fig. [I] shows 
the outage probability as a function of the SNR T, computed 
at each SNR point by evaluating (13i at z = T _1 . Three 
cases were considered for the SINR threshold: (3 — — 10 dB, 




Fig. 1. Conditional outage probability as a function of SNR T. Analytical 
curves are solid, while dots represent simulated values. Top curve: f3 = 10 dB. 
Middle curve: f} = dB. Bottom curve: /3 = — 10 dB. The network geometry 
is shown in the inset. The receiving mobile is represented by the five-pointed 
star at the center of the network, the desired transmitting mobile by the six- 
pointed star immediately above the transmitter, and the 50 interferers are 
shown as dots. 

(3 = dB and (3 = 10 dB. Also shown are results gener- 
ated by simulation, which involved randomly generating the 
exponentially-distributed {gi}. The analytical and simulation 



results coincide, which is what is to be expected because ( 13 1 



is exact. Any discrepancy between the curves can be attributed 
to the finite number of Monte Carlo trials (one million trials 
were executed per SNR point). 

IV. Outage of a BPP 

Because it is conditioned on fi, the outage probability en 
presented in the last section depends on the geometry of the 
particular network, i.e., the location of the interferers. The 
conditioning on ft can be removed by averaging Fz M {z\tt) 
over the spatial distribution of the network. In a BPR a 
fixed number M of mobiles are independently and uniformly 
distributed over the network. Let cm be the spatially averaged 
outage probability when the interferers are drawn from a BPR 
which is found by taking the expectation of Fz M {z\tt) with 
respect to ft: 



e M = E n [e n ] = E [F Zm (T" 1 |n)] = F Zm (T~ 



(14) 



In the above equation, Fz M (z) is the cdf of Zm averaged 
over the spatial distribution, which can be found analytically 
or computed through Monte Carlo simulation. To compute it 
via Monte Carlo simulation, generate a large number N of 
networks, each containing M interferers drawn from a BPP. 
Compute the Fz M (z\H) for each network by using the method 
outlined in Section III, and average over the N networks. 
Letting i~2„ be normalized inverse power coefficients of the 
n th randomly generated network, the Monte Carlo estimate of 



the cdf is 



1 N 

FzA*) = ^E Fz m(*I"»)- (15) 



Note that the Monte Carlo simulation only requires the re- 
alization of the interferer locations, and does not require the 
realization of the fading coefficients. 

For a BPP constrained to an annular region with inner 
radius r ex and outer radius r ne u the spatial coordinates can be 
represented as the complex value Xj = rie^ 01 . The location 
can then be realized by drawing two independent numbers 

x\,i and X2,i from the uniform distribution over {^^j ■> 1 

and [0, 1] respectively and then setting fj = r net ^jx~ij and 
9i = 2irx2,i- 

To avoid the computational burden of a Monte Carlo 
simulation, a closed-form expression for the spatially aver- 
aged outage probability is preferred. Since \\Xi\\ = rj and 
fij = ||A^|| Q , it follows that the normalized inverse power of 
the i th interferer is 



a = 



The pdf of Q,i is 



2 2^o 
— U) <* 

a 



(16) 



(17) 



for r" x < uj < r" et , and zero otherwise. The spatially averaged 
outage probability can now be obtained using ( [D) and ( fl7| ) 
as follows: 

Fz M (z) = J fnHF ZM (z\u>)du (18) 
where the Af-fold integral is over the joint pdf of 



{Qi,...,Qm}- Substituting (13 1 and (17i into (18 1 and using 
the fact that the {f2i} are independent yields: 



Fz M (z) = l 



M 

n 

i=i 

■ 2-Q 



{i- Pi )f3n Q + uj 
(3n + uj 



duj (19) 



Evaluating the integral and assuming pi = p for all users 
results in: 

F Zm (z) = 1- e-^n M K et ) - * K X )} M (20) 
where /3q = /3f2o, 

k = ( r Lt - r lx) 1 (21) 



*(ar) = x" ■ (1 - p) + 



2 • p a;" 



1. 



2 

Q - 



2a + 2 



^ (22) 



and 2-Fi([a, &]; c, x] is the Gauss hypergeometric function 
given by lfT2l 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 





L = 12^-«—" 1 ' *" 


- 


L' = 20 r »*-* ,, ^ ,- *' 




L' = 50^-*-*^*""^ 











30 
M 



Fig. 2. Outage probability ejif as a function of M for five values L' for 
networks with M interfering mobiles drawn from a BPP. The SINR threshold 
is /3 = 3.7 dB, the SNR is set to T = 10 dB, and the other parameters are 
identical to those used to generate Fig. [TJ Analytical curves are solid, while 
dots represent simulated values. 



iFi([a,b];c;x) = 



r(c) 



r(6)r( c - 
i 

b-li 



v 



(l-v) c - b - 1 (l-vx)- a dv (23) 



ii 



where 



T{z) 



t z - x e- l dt. 



(24) 



While the results in this section are for an annular network 
centered upon the reference receiver, we note that other 
network shapes can be accommodated by determining the 
appropriate pdf of the fli and substituting into ( p~8] >. Further- 
more, shadowing can be accommodated by using appropriately 
defined fn.(uj). 

Example #2: Reconsider Example #1, but now instead of 
the network assuming the specific topology shown in Fig. 
[T] let the M interferers be placed according to a BPP. By 
using ( |20| , the spatially averaged outage probability can be 
found. Fig. [2] shows the outage probability as a function of 
M for five values of L' when the SINR threshold is set to 
/3 = 3.7 dB and the SNR is set to T = 10 dB. The values 
of r ex , r net , and a are the same as in Example #1. The 
solid curves show the spatially averaged outage probability 



evaluated analytically, i.e., by using (20 1, while the dots show 
the probability found by Monte Carlo averaging, i.e, by using 
with iV = 10 000 randomly generated networks. Because 
is exact and N large, the analytical and simulation results 



coincide. From Fig. [2] it is observed that the outage probability 
degrades with increasing M and decreasing V . 

Example #3: The dependence on f3 and T is investigated 
in Fig. [3] The values of r ex , r net , and a are the same as in 
Example #1. The number of interferers is set to M = 50 and 
the equivalent number of frequency channels set to V = 200. 
As in Example #2, the spatially averaged outage probability 
is computed analytically and through simulation. In particular, 
the solid curves show the spatially averaged outage probability 
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Fig. 3. Outage probability ej\j as a function of /3 for three values of V 
for an annular network area with outer radius of r ne t = 2, inner radius of 
Train = 0.25. A fixed number of interfering mobiles (M = 50) is drawn 
from a BPP and a path-loss exponent a = 3 is fixed. Analytical curves are 
solid, while dots represent simulated values. 



evaluated using (20 1, while the dots show the probability 
found by Monte Carlo averaging with N = 10 000 randomly 
generated networks. The outage probability is shown as a 
function of the SINR threshold /3 for three values of T. Again 
the analytical curves coincide with the simulation results. From 
[3] it is observed that the outage probability increases with 
increasing j3 and decreases with increasing T. 

V. Outage of a PPP 

Suppose that the network now has a variable number 
of interferers M. Let pm{iti) indicate the probability mass 
function (pmf) of M. Let Fz m {z) be the cdf of Z m when 
there are to interferers drawn from a BPP. It follows that the 
spatially averaged cdf for the variable-sized network is 



Fz(z) = J2p M (m)F Zm (z) 



(25) 



m=0 



and the outage probability averaged over a spatial distribution 
with a variable number of interferers is e = E[€m] = 
Fz (r _1 ), where the expectation is with respect to the distribu- 



tion of M as given by (25 i. We note that the {F Zm {z)} in (25 1 
may be obtained either through simulation using ( p"5| ), which 
requires the summation in ( 25 i to be truncated. Alternatively, 



the analytical expression given by (20i may be used, which 
as will be shown below, does not require truncation of the 
summation. 

When the spatial distribution is a PPP, the distribution of 
M is Poisson with intensity rj = XA, where A is the density 
of the points per unit area and A is the area over which the 
points are distributed. For a PPP of density A, the number of 
interfering mobiles M within area A has pmf 

PM {m) = {XA ^e-* A (26) 



for to > 0. 



Substituting pi and ([26} into ([25l yields 



Fz(z) 



-XA 



E 



(\Ay 



{1 - e-^ {K . [* Ket ) - * (r«J]} m } (27) 



where n and ty(x) are given by (21i and (22 1, respectively. 
By using the identity 1PT31 



Z — ✓ rtl \ 



c • e 



a ■ b 



(28) 



m=0 



and A — it (r^ et — r^j, (27i may be expressed as 

F z (z) = l-cxp{-(3 z-n\(rl et -rl x ) 

x(l -*(»£)])}. (29) 

The expression given in ( |29| generalizes an earlier expres- 
sion given in Baccelli et al. |4| for the outage probability of an 
infinite network (r net — > oo) with no exclusion zone (r ex = 0) 
and constantly transmitting mobiles (p = 1). To see this, set 
p = 1 and r ex = in ((29]) and take the limit as r net —> oo, 

net 



F z {z) — lim 1 — exp < — f3 z — n\r net 



a + 2 



2a + 2 



■/ 1 e t 

AT 



(a + 2) 
By using the identity iTPZl 

2 F 1 ([a, b] ; c, z) = (1 - zy b 2 F 1 [ [b,c - a] ; c, 



z- 1 



(30) 



(31) 



and performing some algebraic manipulations, ( 30 1 becomes 

' ' r a+" 

net 



F z (z) = 1— lim exp < —/3oz 
'a + 2 a + 2 



27rAr Q+2 



2F1 



Because -g 21 

Pa 

simplified to 



/3 (a + 2) V A) 
2a + 2 



1 



00 




1 



(32) 



Fz{z) = 1 — lim exp ^ —f3 z — - — ^^Pq 



1F1 



Since lim 



a+2 a+2 
a a 

r% et /Po 



when r net -+ 00, (|32} can be 

2ttA 

(a + 2)' 

2a + 2 



+ 



1 



(33) 



r„e t ^oo r%JP + 1 

F z {z) = 1 - exp <^ -f3 z - 



1, (33 1 is equal to 

a! 



2 F 1 



2a + 2 



2ir\ 

(« + 2) 
"a + 2 a + 2 
a ' a 
By using the identity lfl2l 

T(c)T(c-a-b) 



a 



-a 



^(Mhd) 



(34) 



(35) 



T(c-a)T(c-b) 
and performing a few algebraic manipulations, d34} becomes 



FzC ) = 1 --xp { ~Pqz - —PIT (-) I fl 2 
' a \ a I V a 



(36) 




Fig. 4. Spatially averaged outage probability e as a function of the mobile 
density A when interferers are drawn from a PPP. Analytical curves are solid, 
while dots represent simulated values. Top solid = 10. Bottom 

solid curve: r ne t = 2. The dotted line curve is the e as a function of A for 
an infinite network. 



which, in the absence of noise, coincides with equation (3.4) 
of JU and equation (61) of OH. 

Example #4: As with Example #2, suppose that ||Xo|| = 1, 
a = 3, and T = 10 dB. Let r ex = and V = 1. The 
interfering mobiles are now placed in a circular region of 
radius r net according to a PPP with node density A. Three 
network radii are considered: r net — {2, 10, oo}. The SINR 
threshold is set to j3 = 3.7 dB. Fig. |4] shows the spatially 
averaged outage probabilities for each of the three values of 
r net as a function of A. The outage probabilities of the two 
networks with finite radius were computed using ( f29] >, while 
the outage probability of the infinite network was computed 



using (36 1. In addition, simulation results are shown for the 
two finite networks, which coincide with the theoretical curves. 
Again, one million trials were executed per A point, and any 
discrepancy between the theoretical result and the simulation is 
due to the finite number of trials. Fig. |4] shows that the outage 
probability increases with increasing A and/or increasing r net . 

Example #5: Fig. [5] investigates the influence of the path- 
loss exponent a and the number of equivalent hopping chan- 
nels L' . As in Example #4, ||X || = 1 and Y = 10 dB. 
The interferering mobiles are placed in an annular region with 
an inner radius r ex = 0.25 and an outer radius r net = 2 
according to a PPP with node density A, and the SINR 
threshold is fixed to f3 = 3.7 dB. Fig. [5] shows the outage 
probability as a function of A for three values of L 1 and three 
values of path-loss exponent a. For each set of (L', a), the 
outage probability was computed analytically using (29 1, as 
shown by the curves, and by Monte Carlo simulation with 
one million trials, as indicated by the dots. Consistent with 
observations made in Examples #2 and #4, Fig. [5] shows 
that the outage probability increases with decreasing L' and 
increases with increasing node density. In addition, Fig. [5] 
shows that the outage probability decreases with increasing 



Fig. 5. Spatially averaged outage probability e as a function of node 
density A for three values U when interferers are drawn from a PPP. The 
source transmitter is placed at unit distance from the reference receiver and a 
exclusion-zone is imposed at the receiver with radius r ex = 0.25. The SINR 
threshold is fixed to f) = 3.7 dB and the SNR is set to T = 10. Analytical 
curves are solid, while dots represent simulated values. 

path-loss exponent a. 

VI. Transmission Capacity 

Often, networks are constrained to ensure that the outage 
probability e does not exceed a maximum outage probability 
C G [0, 1]; i.e., e < £. Under such a constraint, the maximum 
density of transmissions is of interest, which is quantified by 
the transmission capacity (TC) fl4l . With outage constraint (, 
the TC is 

r c (C) = e-^Oa-C) (37) 

where e _1 (C) is the density of the underlying process (BPP 
or PPP) whose spatially averaged outage probability satisfies 
the constraint e < ( with equalitjj^j and (1 — £) ensures that 
only successful transmissions are counted. The TC represents 
the spatial spectral efficiency; i.e. the rate of successful data 
transmission per unit area. With appropriately normalized 
variables, the TC can assume units of bits-per-second per Hz 
per m 2 (bps/Hz/™ 2 ). 

Closed-form expressions for TC can be found in Rayleigh 
fading for networks drawn from either a BPP or a PPP. For 
the BPP case, e _1 (C) is found be solving e = Fz^^r^ 1 ) = ( 



for A. By substituting M = XA into ( |20| i, 

C = l-e-^[ K {*(C t )-*(r«)} 
By solving for A and setting the result to e _1 (^), 

-x (n loga-o + lor- 1 

^ Alog{ K [*(r« et )-*(r«)]}- 
By substituting ^9} into <|37j, the TC for a BPP is 

(l-Q [logq-Q + A.r- 1 ] 

Alog{K[H> (r« et )-f (rgj]}- 



A. I 



(0 



(38) 



(39) 



(40) 



2 Since e is a monotonically increasing function of A, the TC is maximized 
when the constraint e < f is met with equality. 




Fig. 6. Transmission capacity t c (£) as a function of the outage constraint f 
for different values of SNR T, when interferers are drawn from a BPP. The 
network dimensions are r ex = and r ne t = 2. The curves were produced 
using parameters (5 = — 10 dB, V = 1, and a = 3. Top curve: T = 10 dB. 
Middle curve: r = dB. Bottom curve: T = -10 dB. 

Example #6: In this example, a circular network is assumed 
with r ex = 0, r net — 2, and interferers drawn from a BPP. 
The path-loss exponent is a — 3, and the equivalent number 
of frequency channels is L' = 1, which results in p = 1. 
The SINR threshold is set to (3 = -10 dB. Fig. [6] shows the 
transmission capacity as a function of the outage constraint 
£ for three values of SNR V. The curves were produced by 



Fig. 7. Transmission capacity r c (f) as a function of outage constraint £ 
for three network radii, when interferers are drawn from a PPP and r ex = 0. 
The curves were produced using parameters f3 = — 10 dB , V = 10 dB and 
V = 1. Top curve: r ne ± — 2. Middle curve: r ne t — 10. Bottom curve: 

rnet = CO. 

SINR threshold is set to = -10 dB, and the SNR is T = 10 
dB. Fig. [7] shows the transmission capacity as a function of 
the outage constraint £ for the three values of network radius 
r net . The transmission capacity of the two networks with 



using (40 1 and show that transmission capacity increases with 

r. 

When the interferers are drawn from a PPP, e _1 (C) is found 



finite radius was computed using ( 42 1, while the transmission 
capacity of the infinite network was computed using ( |43) . The 
curves show that the TC decreases with increasing network 
radius. 

VII. Modulation-Constrained TC 



by using solving e 
by d29l, 



Fz(r ) = £ for A, where Fz(z) is given 

C = 1 - exp {-^r- 1 - ttXk- 1 {1 - k [* K et ) 
Solving for A and setting the result to e _1 (C), 

io g (i - cr 1 - Por- 1 



L (C) = 



...c- 1 {!-«[* Ket)-* (rgc)]} 
By substituting ((41} into ([37j, the TC for a PPP is 

(1-0 [loga-cr 1 -/^ 

Tc(C) 



TTK 



'- 1 {l-«[*(rS e t)-*(rS B )]} 
When r net —} oo, r ex — and p = 1, d42| becomes 



(0 



(i-o bgtwr 1 -/^ 



-i-i 



This expression agrees with equation (4.10) in 1151 . which 
traces back to equation (62) of lfT4l in the absence of noise. 

Example #7 : In this example, a circular network is assumed 
with r ex = and interferers drawn from a PPP. Three 
different values of r net are considered: r net = {2, 10, oo}. 
The path-loss exponent is a — 3, and the equivalent number 
of frequency channels is L' = 1, which results in p = 1. The 



The transmission capacity expressions presented in the 
previous section are functions of the SINR threshold (3 and 
make no assumptions about the existence of any particular 
* ( r ez)]}} -tyP e of modulation or channel coding. In practice, the SINR 
threshold is a function of the modulation and coding that is 
used. Let C(j) be the maximum achievable rate that can be 
supported by the chosen modulation at an instantaneous SINR 
of 7. If a rate R code is used, then an outage will occur when 
C(7) < R- Since C(j) is monotonic, it follows that is 
the value for which C(/3) = R, and therefore we can write 
P = C-\R). 

Frequency-hopping systems often use noncoherent CPFSK 
modulation (9), ifTTl . The maximum achievable rate of nonco- 
herent CPFSK is given in [9 1 for various modulation indices h, 
where it is called the symmetric information rate. In particular, 
Fig. 1 of J3 shows the symmetric information rate of binary 
CPFSK as a function of 7 for various h. To emphasize the 
dependence of the capacity on h, we use C(h, 7) in the sequel 
to denote the rate of CPFSK with modulation index h. For 
any value of h, the value of the SINR threshold /3 can be 
found from the corresponding curve by finding the value of 
7 for which C(/i, 7) = R. For instance, when R = 1/2 
and h = 1, the required f3 = 3.7 dB. In |16|, it was found 
that in practice, and over a wide range of code rates, turbo- 
coded noncoherent CPFSK is consistently about 1 dB away 



(41) 



(42) 



(43) 



from the corresponding modulation-constrained capacity limit. 
Thus, the /3 required in practice will generally be higher than 
the value of j3 m i n (R, h) by a small margin. For instance, if a 
1 dB margin is used, then the SINR threshold for noncoherent 
binary CPFSK with R = 1/2 and h = 1 should be set to 
/3 = 4.7 dB. 

When accounting for modulation and coding, the maximum 
data transmission rate is determined by the bandwidth B /L of 
a frequency channel, the spectral efficiency of the modulation, 
and the code rate. Let r) be the spectral efficiency of the 
modulation, given in symbols per second per Hz, and defined 
by the symbol rate divided by the 99 percent-power bandwidth 
of the modulation^] The spectral efficiency of CPFSK can 
be found by numerically integrating the normalized power- 
spectral densities given in [11 1, or since we assume many 
symbols per hop, by Equation (3.4-61) of ifTTI and then 
inverting the result. To emphasize the dependence of r\ on 
h, we denote the spectral efficiency of CPFSK as rj(h) in 
the sequel. When combined with a rate-i? code, the spectral 
efficiency of CPFSK becomes Rq(h) (information) bits per 
second per Hz, where R is the ratio of information bits 
to code symbols. The data rate supported by the channel 
is Ri](h)B/L bits per second. The average data rate, or 
throughput, must account for the duty factor d and only count 
correct transmissions. Hence, the throughput is 

Rrj(h)Bd(l - e) Rrj(h)B(l - e 



T 



(44) 



L L> 
The modulation-constrained transmission capacity is the 
throughput multiplied by the node density, 

\Rr)(h)B(l - e) 



r(A) 



XT 



(45) 



In contrast with (37) , this form of transmission capacity 
explicitly takes into account the code rate R, as well as 
the spectral efficiency of the modulation r)(h). It furthermore 
accounts for the hopping bandwidth B/L'. Rather than con- 
straining outage probability, r(A) fixes the node density A 
and allows the outage probability to vary accordingly. Since 



it accounts for the actual system bandwidth B, (45 » assumes 
units of bps/m 2 . By dividing by bandwidth, the normalized 
modulation-constrained transmission capacity 

r'(A) = r_ "WP-«> . m 

takes on units of bps/Hz/m 2 . However, unlike (37i, t'(A) is 



in terms of information bits rather than channel bits. 

VIII. Network Optimization 
by an Exhaustive search 

The main goal of this paper is to find the (L',R,h) 
that maximizes the normalized TC r'(A) for a frequency- 
hopping ad hoc network, assuming that transmissions occur 
using a capacity-approaching code (e.g., turbo or LDPC) and 
noncoherent binary continuous-phase frequency shift keying 

'Percent-power bandwidths other than 99 can be used, but will influence 
the amount of adjacent-channel interference. 



(CPFSK) modulation. The optimization can be accomplished 
using an exhaustive search by performing the following steps: 

1) Pick a value of j3. 

2) Pick a value of h, and determine the rate R corre- 
sponding to the current (3 (this is found by setting 
R = C(h,j3)) and its corresponding bandwidth effi- 
ciency r)(h). 
Pick a value of V . 



3) 
4) 



Use ( 20 1 to compute the average outage probability €m 
if the interferers are drawn from a BPP or ( |25| l to 
compute e if the interferers are drawn from a PPP . 
5) For the set of (h,R) found in step [2] determine r'(A) 
by using ( |46| . 

I 

2 



6) Return to step 

7) Return to step 

8) Return to step 



until all L' are considered, 
until all h are considered. 
1 until all are considered. 
The above procedure will find the t'(A) for each (£', R, h) 
considered, and the optimal value of these parameters are the 
ones that maximize r'(A). By limiting L' to be integer valued 
(which is not necessary if d is a fraction), the number of values 
is finite and an exhaustive search up to some maximum value 
is feasible. The value of (3 is continuous, and therefore must 
be quantized. For the exhaustive search results presented in 
this section, (3 was quantized to a spacing of 0.1 dB over the 
range —2 dB < f3 < 12 dB, and h was quantized to a spacing 
of 0.01 over the range < h < 1. 

A. Optimization Results for a BPP 

The optimization was run for a network of M = 50 
interferers placed according to a BPP with an inner radius 



of r„ 



0.25 and an outer radius r r , 



2. The path-loss 



exponent was fixed to a = 3. Fig. [8] shows the maximum 
normalized modulation-constrained TC T' opt (X) as function the 
SNR T. For each value of T, the optimal set of (L', R, h) that 
maximizes the TC was found using the previously described 
exhaustive search. 

The value of T' opt (\) was computed assuming a capacity- 
achieving code. Suppose that instead, the code has a gap 
of 1 dB from capacity, i.e. that the required threshold /3 is 
1 dB higher than that predicted by information theory. The 
transmission capacity will be lower due to this gap. The curve 
labeled t[ shows the TC of a code when using a code with 
a 1 dB gap from capacity when using the optimal values of 
(L', R, h) found assuming a capacity-achieving code. As can 
be seen, a modest loss in TC occurs when the code has a 
1 dB gap from capacity. In addition, Fig. [8] also shows the 
normalized TC T~' sub (\) of a system with a suboptimal but 
typical choice of parameters: (L',R,h) = (200,1/2,1). The 
results shown in Fig. [8] highlight the importance of parameter 
optimization. The TC is improved by a factor of 5-10 by 
selecting optimal, rather than arbitrary parameters. 



Fig. |9pT| explore the relative importance of each of the three 
parameters. In each curve, the SNR was set to T = 10 dB and 
one parameter is varied. At each value of the parameter, the TC 
is maximized with respect to the other two parameters. Three 
values of path-loss exponent are considered, a — {3,3.5,4}. 
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Fig. 8. Maximum normalized modulation-constrained TC t'(A) as a function 
of the SNR r for three cases: (1) the code achieves capacity; (2) the code has 
a 1 dB gap from capacity; and (3) a typical set of ((£', R, h))=(200, 1/2, 1) 
is used. The interferers are drawn from a BPP. 




Fig. 9. Maximum normalized modulation-constrained TC T^ pt (X) as a 
function of the equivalent number of frequency channels V . The interferers 
are drawn from a BPP and the network dimensions are r ex = 0.25 and r ne t 
= 2. For each value of L', the optimal R and h are found. Top curve: a = 
4. Middle curve: a = 3.5. Bottom curve: a = 3. 



The optimal values of each of the parameters can be identified 
by locating the peaks of each curve. A general trend is that the 
TC improves with increasing a, though the optimal parameter 
values are not strongly influenced by the a. 

B. Optimization Results for a PPP 

Next, the optimization was run for a network with interferers 
drawn from a PPP with an inner radius of r ex = 0.25 and 
an outer radius r net — 2. The SNR was set to T = 10 
dB and the path-loss exponent to a = 3. Fig. 12 shows the 



Fig. 10. Maximum normalized modulation-constrained TC "r' opt {X) as a 
function of the code rate R. The interferers are drawn from a BPP and the 
network dimensions are r cx = 0.25 and r ne t = 2. For each value of R, the 
optimal L' and h are found. Top curve: a = 4. Middle curve: a = 3.5. Bottom 
curve: a = 3. 
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maximum modulation-constrained normalized TC T^ pt (X) as 
function of the mobile density A. For per each value of A, the 



0.004 - 



Fig. 11. Maximum normalized modulation-constrained TC T~' opt (\) as a 
function of the modulation index h. The interferers are drawn from a BPP 
and the network dimensions are r ex = 0.25 and r ne t = 2. For each value of 
h, the optimal L' and _R are found. Top curve: a = 4. Middle curve: a = 3.5. 
Bottom curve: a = 3. 

optimal set of (L',R,h) that maximizes the TC was found 
using the previously described exhaustive search. Similar to 
Fig. [8] Fig. 12 shows the performance r{(A) when the code 
has a 1 dB gap from capacity, and shows the performance 
T 'subW °f a system that uses the typical choice of parameters: 
{V, R, h) = (200, 1/2, 1). While the loss due to using a code 
with a 1 dB gap from capacity is quite minimal, the loss due 
to using suboptimal parameters is quite high, especially in 
sparser networks. 

Fig. [I3|l5l explore the relative importance of each of the 
three parameters by varying one parameter. At each value of 
the parameter, the modulation-constrained TC is maximized 
with respect to the other two parameters. Four values of mobile 
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Fig. 12. Maximum normalized modulation-constrained TC t'(A) as func- 
tion of mobile density A for three cases: (1) the code achieves capacity; 
(2) the code has a 1 dB gap from capacity; and (3) a typical set of 
((£', R, /i))=(200, 1/2, 1) is used. The interferers are drawn from a PPP. 



Fig. 14. Maximum normalized modulation-constrained TC t (X) as a 
function of the code rate R. The interferers are drawn from a PPP and the 
network dimensions are r ex = 0.25 and r n et = 2. For each value of R, the 
optimal V and h are found. Curves from top to bottom: (1) A = 5; (2) A = 
2; (3) A = 0.5; (4) A = 0.1. (4) A = 0.1. 
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Fig. 13. Maximum normalized modulation-constrained TC T„ pt (X) as a 
function of the equivalent number of frequency channels L' . The interferers 
are drawn from a PPP and the network dimensions are r ex = 0.25 and r ne t 
= 2. For each value of L', the optimal R and h are found. Curves from top 
to bottom: (1) A = 5; (2) A = 2; (3) A = 0.5; (4) A = 0.1. 



density are considered, A = {0.1,0.5,2,5}. The optimal 
values of each of the parameters can be identified by locating 
the peaks of each curve. 
In Fig. 



0.65 

h 

Fig. 15. Maximum normalized modulation-constrained TC T~' opt (X) as a 
function of the modulation index h. The interferers are drawn from a PPP 
and the network dimensions are r ex = 0.25 and r ne t = 2. For each value of 
h, the optimal L' and R are found. Curves from top to bottom: (1) A = 5; 
(2) A = 2; (3) A = 0.5; (4) A = 0.1. 



dependence on R, with denser networks requiring a slightly 



13 



it is shown the maximum normalized TC as a 
function of the equivalent number of frequency channels L' for 
a certain annular network area that has three different density 
of interferers per unit Area (A = 5, A = 2 and A = 0.5 ) when 
for each of them the optimal modulation index and the optimal 
rate of the coded modulation are used. Fig. 13 shows a strong 



dependence on L' , with the optimal V becoming larger with 
increasing network density. Fig. 14 shows a relatively weak 



lower rate. Fig. 15 shows a weak dependence on h, with h « 
0.59 providing optimal performance for every density. 

IX. Network Optimization by Gradient search 
method 

The results of the exhaustive search presented in the pre- 
vious section suggest that the modulation-constrained TC is a 
concave function of (L', R, h). It follows that the optimization 
is a convex optimization problem and can be efficiently 
solved through a gradient-search method ifTOl . In particular, the 
optimization can be accomplished by performing the following 
steps: 



TABLE I 

Results of the Optimization for an annular network area 
where the interferers are drawn from a bpp. the number of 
interferers is fixed to m = 50. 





Tex 


a 


T opt 


opi v 




1 


0.25 


3 


0.04427 


0.04427 


66 






3.5 


0.04395 


0.04395 


66 






4 


0.04372 


0.04372 


66 




0.5 


3 


0.04503 


0.04503 


87 






3.5 


0.04468 


0.04468 


66 






4 


0.04437 


0.04437 


66 


2 


0.25 


3 


0.01590 


0.01590 


54 






3.5 


0.01688 


0.01688 


54 






4 


0.01792 


0.01792 


54 




0.5 


3 


0.01641 


0.01641 


63 






3.5 


0.01752 


0.01752 


63 






4 


0.01871 


0.01871 


57 


4 


0.25 


3 


0.00983 


0.00983 


50 






3.5 


0.01187 


0.01187 


47 






4 


0.01395 


0.01395 


50 




0.5 


3 


0.01024 


0.01024 


44 






3.5 


0.01252 


0.01252 


44 






4 


0.01484 


0.01484 


47 



TABLE II 

Results of the Optimization for an annular network area 
where the interferers are drawn from a ppr the intensity a 
per unit area is fixed to a = 1. 



Vnet 


'ex 


a 


' opt 


opt v 


Iv 


l 


0.25 


3 


0.04654 


0.04654 


129 






3.5 


0.04623 


0.04623 


129 






4 


0.04598 


0.04598 


129 




0.5 


3 


0.05932 


0.05932 


129 






3.5 


0.05881 


0.05881 


129 






4 


0.05838 


0.05838 


129 


2 


0.25 


3 


0.01597 


0.01597 


126 






3.5 


0.01697 


0.01697 


126 






4 


0.01801 


0.01801 


134 




0.5 


3 


0.01731 


0.01731 


134 






3.5 


0.01845 


0.01845 


134 






4 


0.01973 


0.01973 


123 


4 


0.25 


3 


0.00977 


0.00977 


117 






3.5 


0.01180 


0.01180 


117 






4 


0.01387 


0.01387 


117 




0.5 


3 


0.01030 


0.01030 


120 






3.5 


0.01258 


0.01258 


135 






4 


0.01491 


0.01491 


123 



1) Pick intervals for L' ([L' min , L' max ]), (3 {{p min , Pmax]) 



' mini hma r 



2) 



3) 
4) 



and h ([h min , h max ])- 
Create sets L set = {L' min ,(L' v 

ftset — {/^mm; ifimax H~~ Pmin/^) i^max\i Slid h se i — 

{hmin, (hmax + h min /2) , h max } composed of the two 
extreme points and center point of each interval. 
Pick one of the three values of /3. 
Pick one of the three values of h and determine the rate 
R corresponding to the current h and (3 (this is found by 
setting R = C(h,j3)) and its corresponding bandwidth 
efficiency rj(h). 
5) For all three values of V and for the set of (h, R) found 



in the last step, determine r'(A) by using (46 1. 
6) Once r'(A) is computed for all three values of L', 
determine which value has the largest normalized TC: 

a) If the maximum is at one of the two external points, 
the center of the search points is moved in that 
direction and two new external points are chosen 
closer to the new center point; 



b) If the maximum is at the center point, the two 
external points are moved closer to the center. Each 
time the maximum is in the center, the distance 
between the two external points and the center is 
gradually decreased. 

7) Return to step [5] and use the new three points, until the 
distance between the two external points and the center 
reaches the values fixed in step[T]and the maximum stays 
in the center. 

8) Repeat step [5] [6] and [7] for all three values of h and save 
the normalized TC of them when L is optimized. 

9) As step [6} once r'(A) is computed for all three values 
of h, determine which value has the largest normalized 
TC: 

a) If the maximum is at one of the two external points, 
the center of the search points is moved in that 
direction and two new external points are chosen 
closer to the new center point; 

b) If the maximum is at the center point, the two 
external points are moved closer to the center. Each 
time the maximum is in the center, the distance 
between the two external points and the center is 
gradually decreased. 

Return to step [8] and use the new three points of h, 
until the distance between the two external points is 
sufficiently small and the maximum point remains in 
the center. 

11) Repeat step [8] [9] and 10 for all three values of f3 and 
save the normalized TC of each. 

12) As step [9) once t'(A) is computed for all three values 
of (3', determine which value has the largest normalized 
TC: 

If the maximum is at one of the two external points, 
the center of the search points is moved in that 
direction and two new external points are chosen 
closer to the new center point; 
If the maximum is at the center point, the two 
external points are moved closer to the center. Each 
time the maximum is in the center, the distance 
between the two external points and the center is 
gradually decreased. 



10) 



a) 



b) 



13) Return to step 11 and use the new three points of /3, 
until the distance between the two external points is 
sufficiently small and the maximum point remains in 
the center. 



The algorithm is initialized by the initial intervals selected 
at step 1. As the algorithm runs, the size of the intervals get 
successively smaller. The algorithm can stop once both of the 
following conditions are satisfied at the same time: (1) the 
optimal values of L', h and /3 are the same of the previous 
iteration; (2) the difference between each of the three elements 
inside the sets L set , h set and j3 se t is equal to the quantization 
step for respectively L, h and f3 . 

Tables [I] and [TT] compare the results of optimizations per- 
formed by exhaustive search and gradient search for networks 



distributed according to BPP and PPP processes, respec- 
tively. The column marked r opt is the maximum modulation- 
constrained TC found by using the exhaustive-search tech- 
nique of Section VIII while T op t v is the value found using the 



gradient-search technique presented in this section. For each 
type of spatial distribution, three values of r net , two values of 
r ex , and three values of a were considered. For the BPP, the 
number of interferers was set to M = 50, while for the PPP, 
the density was set to A = 1. For both processes, the SNR 
was set to T = 10 dB. 

The tables also indicate the number of iterations required 
for the gradient-search technique to converge. The value is 
indicated by the column marked 2y Each iteration requires 
that 200 values of r'(A) be evaluated, since for L was used 
a spacing of 1 over the range 1 < L < 200. Notice the 
slight variation in the number of iterations. This is in contrast 
with the exhaustive-search algorithm, which requires that fixed 
number of values of r'(A) be evaluated. In particular, the 
exhaustive-search optimization considered 2, 848, 200 sets of 
discretized parameters by using the same parameter spacings 
described in Section [Villi 

As we can see from either Table [I] and Table [TT] the gradient- 
search method gives the same maximum normalized TC as 
the exhaustive-search method ( r opt = r optv ). However, the 
gradient-search technique is more efficient because it requires 
fewer values of r'(A) to be evaluated. 

X. Conclusion 

The combination of frequency-hopping, noncoherent 
CPFSK modulation, and capacity-approaching coding is 
a sensible choice for modern ad hoc networks. For such 
systems, the performance depends critically on the number of 
frequency-hopping channels, the modulation index, and the 
code rate. While these parameters are often chosen arbitrarily, 
the system performance can be significantly improved by the 
joint optimization of the three parameters. The modulation- 
constrained transmission capacity is an appropriate objective 
function for the optimization. Preliminary results using 
modulation-constrained TC as the objective function suggest 
that the optimization problem is convex and therefore a good 
candidate for the gradient-search algorithm proposed in this 
paper. 

The derivation of modulation-constrained transmission ca- 
pacity required a careful analysis of the outage probability 
under the assumptions made in this paper. By extending the 
analysis, other fading distributions, such as Nakagami, can be 
considered, as can shadowing. More sophisticated spatial mod- 
els can be considered, for instance by imposing a minimum 
separation among all users. While more sophisticated models 
might not be analytically tractable, they are good candidates 
for the Monte Carlo method proposed in this paper, which 
requires the random placement of mobiles but does not require 
the realization of the fading coefficients. 



The results presented in this paper are just a sample of what 
is possible using this methodology. In addition to considering 
more sophisticated channel models, future work could con- 
sider other network topologies (other than the annular region 
considered in this paper). One example of such a network is 
one where the reference receiver is allowed to move from the 
center of a disk to its perimeter. Other types of modulation 
and reception could be considered, such as nonbinary CPFSK 
with multi-symbol reception [18|. Directional antennas could 
be considered, as could the impact of adjacent-channel inter- 
ference due to the effect of spectral splatter. 
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